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Hydrodynamic interactions that occur between a centrifugal pump impeller 
and a volute are experimentally and theoretically investigated. The theoretical 
analysis considers the inability of the blades to perfectly guide the flow 
through the impeller, and also includes a quasi-one dimensional treatment of the 
flow in the volute. The disturbance at the impeller discharge and the resulting 
forces are determined by the theoretical model. The model is then extended to 
obtain the hydrodynamic force perturbations that are caused by the impeller 
whirling eccentrically in the volute. Under many operating conditions, these 
force perturbations were found to be destablizing. Comparisons are made between 
the theoretical model and the experimental measurements of pressure distribu- 
tions and radial forces on the impeller. The theoretical model yields fairly 
accurate predictions of the radial forces caused by the flow through the 
impeller. However, it was found that the pressure acting on the front shroud of 
the impeller has a substantial effect on the destablizing hydrodynamic forces. 

NOMENCLATURE 


b width of impeller discharge 

h total head (h* = 2h/pfl 2 R 2 ) 

j y^r 

k impeller phase coefficient = cos(tan y ln(R) ) + j sinCtan y ln(R) ) 

r,0 polar coordinate system 

s length in tangential direction 

t time 

v relative velocity in impeller 

w width in volute 

x,y,z rectangular coordinate system 

* Now at Sandia National Laboratories, Albuquerque, NM 87185 
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A.rA.rrA moments of volute cross-sectional area (defined in Equations 14a-e) 
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(i = x,y , j = x,y) t components 2 of generalized hydrodynamic force 
matrix [A] (A„ = A^./pnbQ R 2 ) 

(i = x,y, j=x,yJ components of damping force matrix [C] 

(C? . = C /pnbQR^) 
i J 1 J 

t 2 2 

pressure coefficient at volute inlet = 2(P v (R 2 ,9') -h^/pfl R 2 

* 2 3 

force acting on impeller (F = F/pnbfl R 2 > 
integration constant in Bernoulli's equation 


(i*= x,y , j = x,tf) 
(J ±j = J^/pJtbRj/ 0 ) 

coefficients 

of 

the jerk 

force 

matrix 

[J] 

(i„= x,y , j = x,j), 
(K ij = K ij /pnbn *2* 

components 

of 

stiffness 

force 

matrix 

IK] 

(i*= x,y , j = x,jr) 
(M ij =M i j/p JtbR 2 ) 

components 

of 

inertia 

force 

matrix 

[M] 


* 2 2 

pressure in impeller (P^ = 2P^/pfl R 2 > 

* 2 2 
pressure in volute (P v = 2P v /pQ R 2 > 

impeller radius ( with no subscript, R = R 2 /R 1 > 

radius of pressure tap ring 

velocity in volute ( with no subscript, V* = V 0 ,/fiR 2 ) 
width of impeller at R 2 (W^ = W^/b) 
perturbation function for impeller flow 
angular location of the impeller center (= wt = constant) 
angle of flow path through impeller 

distance between impeller and volute centers (e = e/R 2 > 
fluid density 

2 

flow coefficient = flowrate through pump/2nbJ2R 2 

2 2 

total head rise coefficient = (ha-h^)/pQ R 2 
orbit speed of impeller center (whirl speed) 
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0 rotational speed of impeller (shaft speed) 

Subscripts: 

c,s cos wt and sin wt components (non-dimensionalized) 
d downstream of pump 

exp experimental result 

m force component due to momentum exchange 

p force component due to pressure 

r,0 radial or angular component 

x,y components in rectangular directions (real = x and imaginary = y) 

1,2 impeller inlet and discharge 

Superscripts: 

' measurement made in volute reference frame 

" measurement made from frame fixed to rotating impeller 

* non-dimensionalized quantity 

Special Notation: 

v under bar denotes vector quantity 

V overbar denotes centered impeller value (non-dimensionalized) 

[A] square brackets denote a matrix quantity 

x dot represents a time derivative 

INTRODUCTION 

Several sources, both dynamic and hydrodynamic have been identified as 
contributing to the forces on centrifugal pump impellers. Figure 1 shows a typ- 
ical configuration for a centrifugal pump with a few of the key components 
identified. The primary emphasis of this study was to investigate the forces 
that result from the hydrodynamic interaction between the impeller and the 
volute. The usual design criterion for a volute is that it should provide 
minimum interference to the symmetric impeller discharge flow that would occur 
if no volute was present. However, the discharge flow pattern will depend upon 
the overall flowrate through the impeller. Once the flowrate changes, the 

discharge conditions around the impeller become asymmetric for any given volute. 
Even at the volute design flowrate, the discharge conditions could still become 
asymmetric if the impeller is displaced from the "design" center of the volute 
by shaft deflection, bearing wear, etc.. In either case, the end result of the 
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asymmetric discharge conditions is that there will be a net radial force on the 
impeller (see Figure 1) . 

It is customary in rotordynamic analyses to linearize the radial forces 
acting on the rotor in terms of a steady portion acting on the centered 
impeller, and a time dependent part due to the impeller whirling in a small 
circular orbit. Referring to Figure 1, these terms may be expressed as 



where F and F result from the interaction of the centered impeller with the 
volute, x and th& matrix [A] relates the perturbed force to the eccentric position 
of the impeller. The [A] matrix will be a function of the whirl speed, to, and is 
often expressed as a quadratic in w so that the system resembles a simple stiff- 
ness, damping, and mass model. 

The steady or mean forces, F and F , have been examined in several studies 
and an understanding of them has $een greatly enhanced through papers by Iversen 
et al. [9], Csanady [7], and Agostinelli et al. [2] to name a few. All of these 
authors have shown that there is a particular flowrate where forces on the 
impeller will be minimized for a given volute. Previous experimental [5,10] and 
theoretical [6] investigations have also shown that the components of [A] are 
such that a whirling motion of the impeller would be encouraged rather than 
dissipated by the hydrodynamic effects. This has created concern that the rotor 
assembly may whirl at one of its critical speeds even though the shaft may be 
rotating well above this speed. There also exists the problem of the alternating 
flexural stress that would be developed if the impeller whirls at a 
subsynchronous speed (see Ehrich and Childs [8]). 

In the current study, a theoretical model of the volute and impeller flows 
will be developed and compared to experimental results. Previously, a potential 
flow model for the steady forces on a centered impeller was given by [7] and 
this work was later extended by [6] to include the effects of the impeller whir- 
ling within the volute. Although the potential flow model presents a more clas- 
sical approach of solving for the forces, problems arise in relating the two 
dimensional theoretical volute profile to the three dimensional geometry of a 
real volute. For this reason, a bulk flow description of the flow through the 
volute is chosen for the current work. A similiar treatment of the volute flow 
was presented by [9], but the influence of this flow on the impeller discharge 
conditions was largly ignored and only the non-whirling impeller was considered. 
The impeller/ volute interaction will be included along with the effects of 
impeller whirl in the present analysis. 

THEORETICAL ANALYSIS 

In developing the theoretical model, the problem is broken into its two 
natural parts; models are constructed for the flow through the impeller and in 
the volute. The equations that are generated in these two parts are then 
combined by matching the pressures and velocities at the impeller discharge to 
those at the volute inlet. A full development of this model can be found in 
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reference [1] and only brief summary will be presented here. 

Governing Equations for the Impeller 

Figure 2 illustrates the geometries used in developing the impeller model. 
To relate the pressure between the inlet and discharge of the impeller, a 
simplified unsteady form of Bernoulli's Equation is written as 

P v 2 0 2 r" 2 

r i v 14 r 3v *' 2 r 

~T + ~ 2 ~~ — 2 + J at~ dS “ w e J cos((Dt - at - ©")dr" 

• Q H Q M 

2 r 3 

-we sin(wt - Qt - 0")r"d ©" = F(t) (2) 

s" 

Here the flow is assumed to be two dimensional and the impeller whirl speed 
constant . 

To simplify the model, certain assumptions must be made about the velocity 
field within the impeller. Specifically, the flow in the impeller is assumed to 
follow a spiral path with inclination angle, y, which is fixed relative to the 
impeller for a given flowrate and head rise so that 

= 0" + tan y ln(r"/R 2 ) (3) 

9 9 

Here (r",©") and (R^,©- ) are the coordinates of a general point on a stream- 
line within the imperLer and at the position of discharge respectively. The 
flow path angle, y , of the streamlines is permitted to deviate from the 
impeller blade angle. It is determined in a manner described in Section 2(c) so 
that the theoretical and experimental head/flowrate characteristics coincide. 
To account for the asymmetry caused by the volute, a circumferential perturba- 
tion is superimposed on this impeller flow. This flow perturbation is assumed 
to be stationary in the volute reference frame. Together, these observations 
require that 

v = (v 2 „ + v 2 ,,)^ 2 = 4 Q R 2 ,r " ,Gtt,wt,e)sec y/r" (4) 

The perturbation function p, must from continuity considerations be constant 
along a streamline. For whirl motions with small eccentric orbits, p may be 
linearized as 

p(©",r",£lt,wt,e) =P(© 2 ) + e*{p c (© 2 )cos wt + p s (© 2 )sin wt} (5) 

Equations (4) and (5) can now be substituted into Equation (2). The pressure at 
the impeller discharge is then given as a function of p and the inlet pressure. 
The pressure is not known at the inlet of the impeller, but it can be written in 
terms of the inlet total head which is assumed to be circumferentially constant. 
If there is no pre- swirl at the inlet, this will give the inlet pressure as 
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( 6 ) 


P*(R 1> e i ) = h* - d R p(e 2 ) [d R p(« 2 ) + 2e* J sin(© 2 - wt)} 

- 2 e*d 2 R 2 P(© 2 ) {p c (© 2 )cos wt + p g (© 2 )sin wt} 

for small eccentric displacements. By utilizing Equations (4) through (6) and 
neglecting terms of order e 2 and higher, Bernoulli's Equation can now be 
separated into harmonics with steady, e cos wt, and e sin wt dependence as 

2 dp _ 2 _ 

d sec y[2 ln(R) — + dp ] +D -1=0 (7a) 


2d sec 2 y [ln(R) + d PP C + g ln(R)P s ] + D pc - sin © 2 

+ 2 ^ [d R P sin(© 2 + tan y ln(R) ) - cos(© 2 + tan y ln(R) )/R] 

2 

« 2 
-2 — [cos © 2 - cos(© 2 + tan y ln(R) ) / R] / tan y = 0 (7b) 


dp dD 

U sec rClnlR)^ ♦ * l*f> s - § ln(B)f c ] * D ps + cos e 2 ^ 

- 2 ^ [d R p cos(© 2 + tan y ln(R) ) + sin(© 2 + tan y ln(R) )/R] 

2 

(it) 

-2 — [sin © 2 - sin(© 2 + tan y ln(R) )/R]/tan 2 y = 0 (7c) 


where 


D (©') = D (©') + e* [D (© ' ) cos wt + D (O')sin wt] (8) 

p p pc ps 

In Equations (7a-c) the impeller discharge pressure coefficient, D (O'), has 
been transformed into the impeller reference frame by the approximation, 
©2=0' +e sin(O'-wt). This will prove convenient in the future, because the 
pressure at the impeller discharge is assumed to be equal to that at the volute 
inlet. 


Governing Equations for the Volute 

The geometries used in developing the volute model are shown in Figure 3 . 
The volute flow will be described by a continuity equation, a moment of momentum 
equation, and an equation of motion in the radial direction. Each of these three 
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equations can be written respectively as 


3 (wVq , ) 3(wr'V„,) 

3?' + dr 7 


r ’ ' 

i — = o 


( 9 ) 


3 (wr ' Vq »Vq ,) 
33 ' 


3(wr'r'V 0 ,V r ,) 
3r ' 


+ wr 'r' 


fV 

3t 



p 30' 


( 10 ) 


and 


ap y _ P v e» v e> 
dr' " r' 


(11) 


Here it has been estimated that V and V , (and their gradients) are much less 
than Vq,, except at the inlet of the volute. 

Within the volute, the flow is considered to be primarily in the 0' direc- 
tion and to have a flat velocity profile. This will allow Equations (9), (10), 
and (11) to be integrated over the volute cross-section. When these equations 
are combined with Equations (4) and (5), the pressure and velocity distributions 
in the volute will be given in terms of moments of the volute cross-sectional 
area and the perturbation function, p. Both Equations (9) and (10) can then be 
separated into three parts (steady, e cos a>t, and esinwt) as follows: 

Continuity: 

d(V A) 

jA * - ( 12a) 


d(V A) # d(V cos 0') d(p sin 0 ') 

~lF~ =W i d^ + £sin0' + d[p c + ] (12b) 


d(V A) d(V sin 0 ') 
S — _ w * 

d0 ' i d0' 


0) 

Q 


cos 0' 


+ dtp s 


d(p cos 0 ') 


d0' 


] 


Moment of momentum: 


rA dD d(rAV 2 ) d(V 2 ) 

T~ w = ~ d0^ “ rlnrA + <t(l - dtan y p ) p 


(12c) 


(13a) 
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rA dD. 


DC 


2 d0' 


d(rA V V ) * _ 

- 2 - - 7 - s - ♦ g (W.sin 9-V - rrA V s > 


+ <« + rA) 


d(V^ cos 0') 


d© 1 


- 2rlnrA 


d(V V c ) 
d0' 


+ <tf 2 sin 0 'p 2 + 4 cos 0'(^ + 2 - 24 tan y p)p 

dp W* dD 

+ 4 ( 1-24 tany p)(p„ + sin0' ^) + * 7 "cos 0 ' ^ 


d0 ' 


(13b) 


rA dD PS 

2 d0* 


d(rAVV ) , _ 

" 2 ^ " Q (W i C ° S 9 ' V " ^ V C ) 


~2 


„ d(V sin 0 ') 


+ (VL + rA) 


d0‘ 


- 2rlnrA 


d(V V s ) 
d0’ 


- d 2 cos 0'p 2 + 4 sin 0'(®[+ 2 - 24 tan y p)p 


dP 


dD. 


+ 4 ( 1-24 tan y p) (p a - cos 0' “^p") + ~ sin 


d0' 


(13c) 


where 


P ,3 

A(0 ' ) = j R wdr'/bR 2 , lnrA(0') = J R ln(r '/R 2 )wdr '/bR 2 


rA(0') = J R r'wdr'/bR 2 , rrA(0') = J R r'r'wdr'/bR^ 


rlnrA(0') = j* r'ln(r '/R 2 )wdr '/bR? 


(14a-e) 


and 


V ' 

V*(0') =V(0') + e* [V (0 ') cos wt + V (0')sin wt] = rj" 

C S QR2 


(15) 


In Equations (12a-c) and (13a-c) the perturbation function, p, has been 
transformed into the volute reference frame for convenience in obtaining a solu- 
tion. 


To complete the basic equations for the volute problem. Equation (11) may 
be integrated to give the radial pressure variation in the volute as 
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P*= Dp + 2V 2 [ln(r , /R 2 ) - 8 *cos(©' - wt)] + e*[D pc + 4W C ln(r'/R 2 )]cos »t 
+ e*£D + 4VV ln(r'/R,)]sino)t 

PS S 2 


( 16 ) 


Closure Conditions 

Equations (7), (12), (13), and (16) will describe the flow in the impeller 
and the volute after certain boundary conditions are satisfied. Even though p 
is referred to as the perturbation function, it was never assumed to be small. 
However, from the definition of the flow coefficient, p is required to have an 
average of one. The flow perturbation is further assumed to possess at least 
zeroth order continuity around the periphery of the impeller. This restriction 
on p can be met by satisfying the condition, 

P(R 2 ,0) =p(R 2 ,2n) (17) 


To account for what happens to the volute flow at the tongue, it is assumed 
that the average total head of the recirculated flow will be constant across the 
tongue, that is, 

R 2 (0) Rj(0) 

f, (P + p V~ ,/ 2) I wdr' = f , (P + p V? ,/2) I wdr' (18) 

R 2 (0) ©'= 2n R 2 (0) V * ©'=0 

From the remaining flow that is discharged, the flow path angle, y, will be 
determined. Previously it was stated that this angle will vary with flowrate and 
total head. Using this stipulation, y can be found by equating the predicted and 
experimental total head rises across the pump. This requires that 

iff = = [D~ ( 2n ) + C V 2 ( 2n ) ] / 2 (19) 

exp p v 


where 


C y » 1 + 2[lnrA(2n) - lnrA(O) ] / [A(2n) - A(0) ] 

Admittedly, using an experimental result does limit the preliminary design 
applications of this model. However, the "H/Q" curve (in dimensionless form the 
function ys (c0) is normally available for any pump and it is important that 
this fundamln?al characteristic is properly represented in the model. 

This completes the development of the equations necessary to obtain p, D , 
and V. The nine ordinary differential equations of (7), (12), and (13) were 
solved using centered differencing. The initial conditions of p, D , and V were 
chosen in an iterative manner to satisfy the closure conditions stated above. 
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Hydrodynamic Forces on the Impeller 


Basically, there are two sources that contribute to the radial hydrodynamic 
forces on the impeller. One part is due to asymmetric pressure distribution 
around the impeller. The other is caused by the asymmetric momentum fluxes at 
the impeller inlet and discharge. The first contribution is evaluated by 
integrating the pressure around the inlet and discharge of the impeller: 


.2 71 




2n 


0 "l-xJo p l<W R 2 e 




'de„ 


( 20 ) 


where j denotes the imaginary part that corresponds to the y direction (see Fig- 
ure 2) . The second contribution is found by applying the momentum equation to 
obtain 


F (F + jF ) 

~m = x J v m _ 

pb pb 


3 2n ^2 

_ejQt 3t Jo J Ri <v r "+ jv e ")e je "r"dr"dO" 

2tt ^2 

- e jnt [J 0 (v r "+ jv e ")v r "e je ”r''dO"] I 

R 

- e^ t 2Qj* 0 J R (jv r "- v e ")e je "r”dr"dO" 

+ u> 2 e7r(R 2 - Rj)e jwt 


( 21 ) 


When the pressure distributions of Equations (6) and (7) and the velocity 
profiles described by the no inlet pre-swirl condition and Equations (3) — ( 5) are 
applied to Equations (20) and (21), the resulting force on the impeller is 

F =F +F=F + e (F cos wt + F sin cat) (22) 
~ ~p ~m ~ ~c ~s 


where 

2tt J0 

F = rf 2 [W*sec 2 y + k R - 2 + 2 j tan y ] J Q p 2 (« 2 )e 2 dO/27t 

* 2 r 2n - J6 2 

- jdtW.sec y ln(R) + ll J Q p(« 2 )e dO^Tt 


2 * 9 

F q = 4 [W i sec y + kR-2 + 2j tan y] j 0 P(^ 2 )P c (0 2 )e d0 2 /7T 

* 2 f 2lt 

- jdfW^ec y ln(R) + 1] J Q P c (« 2 )e dO^n 

hi * 2 r 2n 

+ q d[W i sec y ln(R) + k/R-l]J 0 P s (^ 2 )e 

w * r 2n - j6 2 

+ 0 ^ ® Jq P(0 2 )sin(0 2 + tan y ln(R))e d0 2 /rr] 


(23 a) 
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(23b) 


w 2 

- q[2j<* + W*/(k R) ] - — (W*[l - 1/ (k R)]/tan 2 r - 1 + 1/R 2 } 

F g = rf 2 [W*sec 2 y + k R - 2 + 2j tan r] 

* 2 r 2n 

- j^[W i sec yln(R) + 1] J Q P s <« 2 )e d^/n 

- ® d[W*sec 2 y ln(R) +k/R-l]J 0 P c (e 2 )e j 2 d0 2 /jr 

u * r 2n - j0 2 

" q WHWJ 0 p(« 2 )cos(« 2 + tany ln(R))e d^/n] 

w 2 

- q j [ 2 jaf + W*/(k R)] - — j {W*[l - l/(k R) ] /tan 2 y -1 +1/R 2 } (23c) 

a 


and, k cos(tan y ln(R) ) + j sin(tan y ln(R) ) . Expressed in the terms used in 
Equation (1), these components are 


F = F + jF 
~ x J y 


~c A xx + ^ A yx 


and F = A* + jA* 
~s xy ° yy 


(24a-c) 


Presentation of the calculated results will be postponed so that the 
experimental and theoretical results can be discussed together. 

TEST FACILITY 


The experimental results presented in this paper were obtained using the 
Rotor Force Test Facility at the California Institute of Technology, Pasadena. 
Details of the equipment have been given in previous papers [3,4,101, so only a 
brief description will be presented here. Figure 4 shows the test section where 
the centrifugal pump being examined is located. The impeller is mounted on the 
internal balance and the entire assembly is turned by the main shaft. The main 
shaft passes through an eccentrically drilled cylinder, which when rotated, 
causes the impeller to whirl in a 0.0990 inch diameter circular orbit. Forces on 
the impeller are sensed through strain gauges on four posts located in the 
internal balance. The relationships between the strains and forces were found by 
static calibration tests. 


Descriptions of the impeller and one of the volutes that were tested are 
given in Figures 5 and 6. The impeller (referred to as Impeller X) is a five 
bladed cast bronze impeller with a specific speed of 0.57 and blade angle of 
65®. The 86® spiral volute (Volute A) is constructed of fiberglass and designed 
to be "well matched" with Impeller X at a flow coefficient of 0.092. The dimen- 
sions of the volute cross-sections, shown in Figure 6, were used in evaluating 
the integrals of Equations (14a-e). 

Two modifications have been made on the test facility for the benefit of 
this research. They were considered necessary in order isolate the interaction 
between the impeller and the volute from external influences. The modifications 
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are illustrated in Figure 7 . To separate the flow in the volute from the annular 
gap region, rings were installed 0.005 inch from the edge of the impeller. In 
addition, the flange of the test section was removed so that the front shroud of 
the impeller was exposed to the "reservoir-like" conditions of the test chamber 
(see Figure 4) . 

The removal of the front flange of the test section was judged to be essen- 
tial after pressure measurements were made in the annual gap region with the 
flange in place and the rings removed. The measurements indicated that the fluid 
trapped in this region was responsible for a hydrodynamic stiffness (see Equa- 
tion (26)) given approximately by 
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When compared with Chamieh's [5] direct measurements of the total hydrodynamic 
stiffness on the impeller (annular gap plus volute) given approximately by 
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it is seen that the contribution from the annular gap is significant. With the 
flange removed, it was anticipated that the fluid forces on the front shroud of 
the impeller would be largely eliminated. 

COMPARISONS BETWEEN EXPERIMENTAL AND THEORETICAL RESULTS 

A preliminary step in the theoretical calculations must be the estimation 
of the impeller flow path angle, y (see Section 2(a)). In practice, information 
on the actual total head rise as a function of flowrate is almost always avail- 
able; an example for Impeller X and Volute A is presented in Fig. 8. By setting 
4j=\ jj , the flow path angle, y, shown in Fig. 9 was obtained. Note that the 
typical magnitude of y is about 80° while the blade angle of Impeller X is 65° . 

Measurements of the static pressure of the discharge from the impeller were 
made using holes drilled at the inlet to the volute (see Figures 6 and 7) . The 
circumferential pressure distributions are compared with the theoretical results 
in Figures 10 and 11. The pressure taps were alternately placed in the front and 
back of the volute, resulting in the slight oscillation of the data. The 
results were obtained for a range of shaft speeds from £2 = 800 to 1200 RPM, but 
the non-dimensional ized pressures were found to be independent of the speed. 
Figure 10 shows that the theory gives a good approximation of the the pressure 
distributions over a moderate range of flow coefficients. For flow coefficients 
larger than this range, the correlation begins to falter as shown in Figure 11. 
It was concluded that the deviation was caused by the inadequacy of a one dimers 
sional treatment of the flow near the tongue of the volute. At the higher 
flowrates, it has been suggested [11] that there is a reversal of the direction 
of flow in the region just inside the tongue. The effect on the pressure 
distribution of displacing the impeller is also demonstrated in Figure 11. The 
model appears to follow the changes that occur, even when the absolute pressure 
predictions are rather poor. 
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A comparison between the experimental and theoretical steady forces on the 
impeller is given in Figure 12. One set of experimental results was obtained by 
placing the impeller in four equally spaced orbit positions and then averaging 
the internal balance force measurements. The second set (for 4-0.06 and 0.10) 
was obtained by integration of the discharge pressure measurements. The theoret- 
ical model tends to overpredict the steady or average radial forces somewhat, 
but it does give reasonable results considering the crudeness of the model. 
Coldlng-Jorgensen's [6] steady force calculations for a 67.5° blade angle 
impeller in an 86° spiral volute are also shown in Figure 12. The present model 
appears to give a more accurate assessment of the steady forces when compared 
with the experimental results. The agreement between the two sets of experimen- 
tal data indicates that the primary cause of the radial force is the asymmetric 
pressure distribution at the discharge of the impeller. Moreover, the theoreti- 
cal model predicted that the discharge pressure was responsible for 99% of the 
total force on the impeller while the net momentum flux contribution was essen- 
tially negligible. It might also be of interest to note that over the entire 
range of flowrates for which theoretical results are presented, the predicted 
perturbation in the impeller discharge flow never exceeded 6% of the mean flow. 

Figure 13 presents the components of the generalized hydrodynamic force 
matrix, [A] , that result when the impeller whirls in an eccentric orbit at the 
pump design flowrate (d = 0.092). From the experimental data, it is seen that 
the cross-coupled terms (i.e. A , A ) imply that forces act in the direction 
of the whirl orbit up to ft>/n = 0.T0. This destabilizing influence is predicted 
by the theoretical model to occur up to o>/G = 0.14. Due to the coupled nature 
of Equations (7b) and (7c), it was not possible to calculate [A] beyond the 
range of whirl ratios shown in Figure 13. This problem is believed to be the 
result of the current limitations of the iterative technique used in obtaining 
the solution. 


As was mentioned in the introduction, it is a standard practice to express 
the matrix elements of [A] in powers of e. By examining the A term in Figure 
13 , it is apparent that a quadratic in w will not adequately describe the 
features of the matrix element. A cubic, however, can approximate all of the [A] 
matrix element variations with w giving the coefficients of such an expansion 
as 
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or ail ter natively as 
_ x 
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x = e cos iot and y = e sin wt 


The [K] , [C] , and [M] matrices correspond to the stiffness, damping, and iner- 
tial components that are commonly employed in rotordynamics. Since the [J] 
matrix is related to the third order time derivative of the impeller displace- 
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ment (which is conventionally known as the jerk) , it will be referred to as the 
"jerk" matrix. 

The resulting [K] matrix elements of the cubic expansion are given in Fig- 
ure 14, and the [C] , [M] , and [J] matrix elements are presented in Figure 15. 
Included in Figure 14 are the stiffnesses that were calculated using the force 
measurements (from the internal balance and the pressure distributions) taken at 
four impeller positions. Also shown in Figure 14 are the stiffnesses predicted 
by Colding-Jorgensen [6] for an 86° spiral volute. With the exception of the K 
term, the current model does a fair job of describing the variation of stiffness 
with flowrate. The magnitudes, however, tend to be underpredicted by the theory. 
Over most of the range of flow coefficients, the stiffness is such that it would 
encourage the whirling motion of the impeller. The same is also true of the 
damping when the flowrate drops below i - 0.07 as shown in Figure 15. The 
magnitude of the damping components computed by [6] (not shown here) were less 
than 10% of those predicted by the present model. In general, the inertial 
force would discourage an orbital motion of the impeller, but it will tend to 
drive the impeller in the direction of the displacement. The jerk force attains 
significant values only at the lower flow coefficients. 

CLOSING COMMENTS 

A theoretical model has been developed to describe the flow in the impeller 
and the volute, along with the interactions that occur between them. This 
investigation was undertaken to provide a better understanding of the 
destabilizing hydrodynamic forces that have been observed [5,10] on a whirling 
centrifugal pump impeller. To implement the model requires only a knowledge of 
the dimensions of the volute and impeller, and the total head rise across the 
entire pump. Comparisons between the predicted and experimental results are 
encouraging. Experimentation with different volute geometries and over a wider 
range of operating conditions (flow coefficient and whirl ratio) would provide a 
more crucial test of the theoretical model. It might also prove insightfull to 
incorporate the effects of inducers and diffuser vanes into the theoretical 
model. These devices are now commonly employed on many high performance 
centrifugal pumps. 

Previous experimental results [5,10] have tended to over-estimate the 
contribution of the volute /impeller interaction to the total stiffness force 
acting on the impeller. The over-estimation came about because of an asymmetric 
pressure distribution in the fluid trapped on the front shroud of the impeller. 
Since real pumps do have fluid in this region, it will be important in the 
future to perform a detailed study of this area. 
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Figure 1. Description of a centrifugal pump. 
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Figure 2. Geometry of a centrifugal pump impeller. 



Figure 3 . 


Geometry of a volute. 










Figure 4. Schematic of the test section. 



Figure 5. Dimensions of Impeller X. 
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Figure 6. Dimensions of Volute A. 
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Figure 7. Modifications of the test facility. 
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Figure 8. 


Figure 9. 



Flow Coefficient, <p 

Total head coefficient vs. the flow coefficient for the Impeller X 
and Volute A pump. The solid line represents the average value of 
the experimental data. 



Flow Coefficient, <£ 

The average angle the flow follows through the pump (measured with 
respect to the radial line) as a function of flowrate. 
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Figure 10. Experimental and theoretical pressure F1 8 ure 11 • Experimental and theoretical pressure 

distributions at the inlet of the volute distributions at the inlet of the volute 

for three different flowrates with the f or two different flowrates with the 

impeller placed near the tongue. Measure- impeller placed near to and far from the 

ments were taken with Q = 1000 RPM. volute tongue. 
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Figure 14. Hydrodynamic stiffnesses as functions of the flowrate. Experimental 
results are from internal balance and pressure measurements. 
Colding-Jorgensen's [6] results are for an 86° spiral volute with a 
67.5° blade angle impeller. 
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Flow Coefficient, <p 


Figure 15. Damping, mass, and jerk force coefficients as functions of the 
flowrate as predicted by the current theoretical model. 
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